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LONG-TERM  GOALS 

The  long-term  goal  of  this  project  is  to  understand  the  role  of  ocean  turbulence  in  submesoscale 
dynamics. 

OBJECTIVES 

The  objective  of  this  project  is  to  examine  the  role  of  ocean  turbulence  in  submesoscale  mixing 
observed  during  the  June  2011  LatMix  experiment.  Key  questions  to  be  addressed  are  the  following: 

(1)  What  are  the  space  and  time  statistics  of  breaking  internal  waves  and  how  often  do  they  result 
in  turbulence? 

(2)  What  is  the  role  of  fine  scale  intrusions  in  along  and  across  isopycnal  mixing?  Is  this  the 
pathway  of  a  forward  cascade  of  energy  and  scalar  variance  from  the  submesoscale? 

(3)  What  is  the  nature  of  the  scalar  spectra  along  isopycnals  (dye,  spice)  between  the 
submesoscale  and  the  microscale  ? 

BACKGROUND 

Central  to  quantifying  ocean  mixing  is  the  relationship  of  essentially  two  dimensional  QG  like 
submesoscale  variability,  which  produces  along  isopycnal  mixing,  to  that  of  three  dimensional 
turbulence,  which  produces  across  isopycnal  mixing.  Strongly  related  to  this  is  the  nature  of  the  scalar 
variance  spectral  wavenumber  cascade.  Numerical  studies  such  as  that  of  Molemaker  et  al.  (2005), 
and  Thomas  et  al.  (2007)  suggest  a  forward  scalar  variance  cascade  within  the  submesoscale  due  to 
mesoscale  straining.  This  is  hypothesized  to  result  in  a  “Batchelor”  like  -1  power  law  for  the  scalar 
variance  wavenumber  spectrum.  Thus,  a  key  factor  in  addressing  the  three  original  LatMIX 
hypotheses  and  in  quantifying  lateral  along  isopycnal  mixing,  is  the  nature  of  the  scalar  wavenumber 
spectra  and  the  scalar  variance  cascade  process.  However,  it  should  be  noted  that  a  very  significant 
assumption  is  being  made  here,  namely  that  there  is  stationarity  and  homogeneity  over  some  (as  of  yet 
undetermined)  space  and  times  scale.  An  alternative  viewpoint  is  that  meso/submesoscale  shear  and 
straining  can  result  in  sporadic  small  scale  intrusions  which  can  lead  to  the  pathway  of  small  scale 
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diffusion  and  dissipation.  Thus,  a  spectra  approach  at  the  submesoscale  must  be  carefully  formulated 
to  take  this  into  account.  In  addition,  care  must  be  taken  to  isolate  the  role  that  breaking  internal  waves 
play  in  this  connection  of  submesoscale  dynamics  to  that  of  turbulence. 

A  remarkable  accomplishment  of  quasi  equilibrium  turbulence  theory  has  been  the  development  of  the 
-5/3  law  and  -  3  law  for  the  kinetic  energy  wavenumber  spectra  of  3D  and  2D  (QG)  turbulence, 
respectively.  These  laws  are  based  on  ideas  originally  developed  by  A.  N.  Kolomogorov  (1941)  and 
then  later  by  G.  K.  Batchelor  (1953).  Their  experimental  and  numerical  verification  is  one  of  the  great 
scientific  accomplishments  of  the  past  several  decades.  These  spectral  formulations  are  based  on  a 
combination  of  nonlinear  cascade  physics  and  dimensional  analysis  (Vallis,  2006).  In  the  3D 
turbulence  case,  the  rate  of  transfer  of  kinetic  energy  (KE)  is  constant  and  equal  to  s,  the  turbulent 
dissipation  rate,  with  non  linear  3  D  processes  resulting  in  KE  being  cascaded  in  the  forward  direction 
to  smaller  scales  and  higher  wavenumbers.  Note,  however,  that  s  has  two  separate  meanings.  It  is  both 
the  amount  of  kinetic  energy  which  is  eventually  dissipated  and  it  is  also  the  rate  of  transfer  of  kinetic 
energy  from  larger  turbulent  scales  to  the  mid  range  turbulent  scales  and  then  on  to  the  smaller  scales 
(Tennekes  and  Lumley,  1972).  For  the  2D  QG  turbulence  case,  as  a  result  of  potential  vorticity 
conservation,  non  linear  processes  result  in  enstrophy  being  cascaded  in  the  forward  direction  at  some 
constant  rate,  \\i.  This  results  in  KE  being  cascaded  in  the  inverse  direction  to  larger  scales  and  lower 
wavenumbers.  Thus,  \\ /  in  QG  turbulence  plays  the  same  type  of  role  as  a  constant  cascade  parameter 
as  s  does  in  3D  turbulence.  It  is  straightforward  to  show  that  the  rate  of  transfer  of  KE  in  QG,  2  D 
turbulence,  s,  is  not  constant  over  this  wavenumber  regime. 


These  same  quasi  equilibrium  cascade  concepts  can  be  applied  to  the  scalar  spectra.  Let  %  be  the  rate  of 
transfer  of  some  scalar  such  as  temperature  and  salinity.  If  we  assume  an  equilibrium  subrange  in 
which  there  is  an  imposed  larger  scale  strain  rate,  S,  on  the  scalar  spectra,  then  %  must  be  constant  in 
this  wavenumber  subrange.  From  just  dimensional  scaling  the  wavenumber  scalar  spectra  O  must  take 
the  form 


(1) 


d>  =  A—k~x 
S 


for 


K  <K  <K, 


(Kraichnan,  1967),  where  A  is  a  universal  constant  of  order  1;  x ,  the  rate  of  transfer  of  scalar 

variance,  has  units  ^sca*ar  un^s)  an(j  g  jiag  unjts  0f  (time)'1  .  In  this  analysis  we  will  take 

time 

spectra  as  referring  to  the  ID  (measurable)  horizontal  wavenumber  spectra. 


our 


In  the  3D  turbulent  case  the  inertial  subrange  strain  rate,  S,  must  scale  with  the  constant  transfer  rate  of 
kinetic  energy,  s,  and  only  local  wavenumber  (Batchelor,  1953,  Tennekes  and  Lumley,  1972),  resulting 
in 

_5 

(2)  d>  =  A—k~3  Ke<K<Kv 

£ 

while  in  the  diffusive  subrange  S  must  scale  with  the  local  strain  rate  produced  by  the  larger  scale 
eddies  in  the  inertial  subrange  resulting  in 
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(3)  =  K,<K<Ka 

V 

where  Ke ,  is  the  energy  containing  wavenumber  limit,  kv  the  Kolomogorv  wavenumber  limit,  and  kb 
,the  Batchelor  wavenumber  limit.  Note  that  the  constant  B  can  be  obtained  in  terms  of  A  by  setting  the 
two  forms  of  the  wavenumber  spectra  given  by  (2)  and  (3)  equal  at  the  wavenumber  boundary  kv  .  The 

well  known  and  well  verified  spectral  forms  (2)  and  (3)  represent  the  scalar  inertial  and  diffusive 
subranges  produced  by  3D  turbulence  (Dillon  and  Caldwell,  1980). 

These  ideas  can  be  extended  to  the  submesoscale  regime  for  the  case  of  variability  modeled  as  2D  QG 
turbulence.  Taking  S  as  the  mesoscale  strain  rate  we  see  that  equation  (1)  then  describes  the  scalar 

spectrum  in  the  submesoscale  which  can  be  thought  of  as  a  “Batchelor”  spectra  with  the  mesoscale 

1 

£  — 

strain  rate  S  replacing  the  3D  turbulent  strain  rate  (— )2  of  equation  (3).  Recent  results  of  Ferrari 

v 

(ONR  LatMIX  Boston  workshop,  2008)  and  others  (Rudnick,  2001)  suggest  that  “spice”  spectra  in  the 
submesoscale  follow  the  -1  power  law  of  equation  (1)  and  support  theoretical  concepts  of  Molemaker 
et  al.  (2005)  and  Thomas  et  al.(2007)  and  that  it  is  the  mesoscale  strain  rate  S  which  results  in  the 
scalar  variance  forward  cascade  into  the  submesoscale. 

Internal  wave  induced  vertical  displacements  of  a  scalar  gradient  such  as  that  of  temperature  and 
salinity  can  also  result  in  scalar  variability.  One  of  the  great  accomplishments  of  the  past  30  years  has 
been  the  near  universal  nature  of  the  Garrett  and  Munk  (1975)  model  in  describing  IW  spectral  forms, 
although  it  should  be  noted  that  very  limited  horizontal  spectral  measurements  of  internal  wave  have 
been  made.  The  GM  models  predict  that  scalar  variance  spectra  would  go  as 

(4)  <J>  =  E^k~2  for  k{}  <k  <  kc  , 


where  k0  ,  kc  are  the  wavenumber  limits  of  the  internal  wave  filed,  determined  by  the  source  and  sink 


of  internal  wave  energy.  In  equation  (4)  E,  has  units  of 


(scalar  units)2 


One  of  the  most  important 


- - - - - OJ  ■  ■■■  - V  V  - -  .  -1  •  - - - -  f 

length 

recent  theoretical  and  observational  results  is  in  relating  GM  type  internal  waves  to  turbulent 
dissipation  rate,  s  (Henyey,  1986,  Gregg,  1989,  Polzin,  1995)  and  that 

(5)  e  oc  E  2 


This  is  related  to  the  fall  off  of  the  spectra  beyond  k  >  kc  where 


E,  . 

(6)  o=^r3 

K 

From  vertical  microstructure  profilers  strong  observational  support  for  (5)  and  (6)  occurs  for  the 
vertical  ID  wavenumber  spectra  of  temperature.  No  conclusive  set  of  observations  have  shown  the 
validity  of  (6)  for  the  horizontal  wavenumber  spectrum. 
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However,  recently  Klymak  and  Mourn  (2007a, b)  did  attempt  to  examine  the  nature  of  the  horizontal 
wavenumber  spectra  over  the  internal  wave  and  submesoscale  ranges  of  1  m  to  1  km.  They  found 
some  limited  evidence  of  (6)  but  for  the  most  part  observed  that  the  horizontal  wavenumber  spectrum 
seemed  to  follow  equation  (2).  This  is  an  unexpected  result  since  in  this  wavenumber  regime  3D 
turbulence  is  expected  to  be  highly  anisotropic,  which  some  have  argued  in  analogy  to  the  atmospheric 
case  (Lumley,  1964,Weinstock,1985,  Holloway,  1986)  should  result  in  a  “buoyant”  subrange  with  the 
spectra  following  a  -3  power  law  as  in  equation  (6).  Moreover  the  type  of  measurements  performed  by 
Klymak  and  Mourn  (2007a,  b),  towing  an  instrumented  body  at  1  m/sec  at  essentially  constant  depth, 
are  susceptible  to  finestructure  contamination.  This  would  result  in  flattening  a  -3  power  law  into  a 
power  law  closer  to  -5/3  (Goodman,  1978). 

Performing  scalar  measurements  such  as  temperature  and  salinity  along  an  isopycnal  (spice)  and  at  the 
same  measuring  the  displacement  of  the  isopycnal  and  the  turbulent  fields  would  be  a  clear  way  of 
distinguishing  these  mechanism  and  sorting  out  the  different  contributions  of  submesoscale  variability, 
internal  waves,  and  3D  anisotropic  turbulence  in  producing  along  and  across  isopycnal  mixing,  which 
is  critical  to  examining  the  nature  of  submesoscale  mixing. 

APPROACH 

The  observational  approach  has  been  to  use  the  Autonomous  Underwater  Vehicle,  T-REMUS,  shown 
in  Fig.  1.  T-REMUS  is  a  custom  designed  REMUS  100  vehicle  manufactured  by  Hydroid  Inc., 
containing  the  Rockland  Microstructure  Measurement  System  (RMMS),  an  upward  and  downward 
looking  1 .2  MHz  ADCP,  a  FASTCAT  Seabird  CTD,  and  a  WET  Labs  BB2F  Combination  Spectral 
Backscattering  Meter/  Chlorophyll  Fluormeter.  For  the  2011  LatMiX  experiment  the  BB2F  system 
was  configured  with  a  Fluorescein  sensor  to  detect  the  fluorescein  dye  deployments.  In  addition,  the 
vehicle  contains  a  variety  of  “hotel”  sensors  which  measure  pitch,  roll,  yaw,  and  other  internal 
dynamical  parameters. 

This  suite  of  sensors  on  T-REMUS  allows  quantification  of  the  key  dynamical  and  kinematical 
turbulent  and  submesoscale  physical  processes.  The  turbulence  measurements  are  made  concomitantly 
with  very  high  spatial  resolution  measurements  of  velocity  shear,  temperature,  salinity,  and  depth.  The 
parameters  which  can  be  estimated  from  the  data  collected  by  the  T-REMUS  include:  the  turbulent 
velocity  dissipation  rate,  s,  the  temperature  variance  dissipation(diffusion)  rate,  x,  the  local  buoyancy 

du 

frequency,  N,  and  vertical  shear,  — ,  and  thus  using  the  latter  two  quantities  the  Richardson  (Froude) 

dz 

N2  _2 

number  Ri  =  — - —  =  Fr  .In  addition  we  obtain  simultaneously  the  local  submesoscale  fields  of 

.du.  2 

dz 

velocity  salinity,  density,  fluorescein,  chlorophyll-a  and  optical  (700nm)  scattering  (Goodman  & 

Wang,  2009;  Wang  &  Goodman,  2009)  in  which  the  turbulence  is  embedded. 
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Figure  1,  The  SMAST  T-REMUS  Autonomous  Underwater  Vehicle.  It  is  2.0  m  long,  20  cm 
diameter,  and  63kg  mass.  Vehicle  based  sensors  are  indicated  in  the  figure. 


WORK  COMPLETED 

We  carried  out  a  very  successful  series  of  deployments  of  T-REMUS  during  the  June  LatMix  20011 
experiment.  In  Figure  2  we  show  a  schematic  of  the  T-REMUS  deployments.  The  T-REMUS  was 
operated  to  make  a  series  of  box  patterns  around  a  drifting  buoy  drogued  at  35  m  depth.  The  idea  was 
to  follow  the  submesoscale  flow  field  while  simultaneously  obtaining  fine  and  microstructure 
(turbulent)  data. 


Figure  2,  Schematic  map  of  the  T-REMUS  deployments  for  the  LATMIX  DRI 2011  experiment.  The 
T-REMUS  AUV  follows  the  Gateway  Buoy  by  means  of  an  attached  drogued  buoy . 
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To  accomplish  this  type  of  deployment  and  track  and  communicate  to  and  from  the  vehicle  from  our 
supporting  ship,  the  Oceanus,  we  developed  a  unique  three  buoy  system,  shown  in  Figure  3.  A  leader 
buoy  with  a  10  m  holey-sock  drogue  set  between  30  to  40  meters  was  used  to  follow  the  flow  field  (be 
Lagrangian)  at  the  approximate  depth  range  of  the  dye  deployments  during  LatMix  2011.  The  second 
buoy  contained  a  thermistor  chain  set  from  1  meter  depth  49  meter  depth  containing  16  equally  spaced 
therimsitors.  An  RBR  CTD  was  set  at  the  end  of  the  T-chain  and  a  depth  sensor  at  the  top  of  the  chain. 
The  third  buoy  contained  the  Hydroid  communication  "gateway"  system.  It  was  set  to  provide  acoustic 
link  from  the  buoy  to  the  vehicle  and  then  RF  to  the  ship.  We  have  subsequently  added  Iridium 
satellite  communication  capability. 


Figures  3  Three  buoy  drogue  Gateway  Communication  system  used  in  LatMix  2011.  The  leader 
buoy  has  a  10  m  holey-sock  drogue  centered  at  10  m  depth.  The  second  buoy  contains  a  16  element 
thermistor  chain  spanning  lm  to  49  m.  The  third  buoy  is  the  hydroid  Gateway  communication  buoy. 


In  Figure  4  we  show  the  times  and  locations  of  the  deployments  made  between  3  June  2011  and  18 
June  2011.  Also  shown  are  the  actual  tracks  of  the  T-REMUS,  both  in  absolute  and  relative 
coordinates. 

Data  Analysis  and  Preliminary  Results 

Much  of  the  work  following  the  2011  field  deployments  was  in  reducing  the  enormous  amount  of  data 
taken  to  a  usable  scientific  format.  As  of  this  writing,  over  50  %  of  this  phase  has  been  completed. 
Preliminary  analysis  has  been  initiated  on  the  second  half  of  the  2011  LatMix  experiment  when  the  T- 
REMUS  via  our  three  buoy  system  intruded  into  a  frontal  zone.  We  have  also  written  a  number  of  new 
analysis  routines.  We  have  developed,  following  Winters  and  D'Asaro  (1996),  a  three  dimensional 
sorting  algorithm  to  obtain  turbulent  patch  statistics.  We  will  present  below  some  of  our  results  to  date. 
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LatMix  T-REMUS  Deployments 


Date  in 

06/03/1 1 

06/05/1 1 

06/08/1 1 

06/10/11  * 

06/14/11  - 

06/16/11  - 

06/18/11 

Time  in  (UTC) 

06:03:14  PM 

02:34:42  PM 

05:46:34  PM 

11:32:09  AM 

01:17:11  PM 

01:03:08  PM 

01:03:24  PM 

Time  out  (UTC) 

08:25:23  PM 

10:50:43  PM 

11:16:32  PM 

09:44:37  PM 

10:57:00  PM 

10:58:05  PM 

08:23:43  PM 

Starting  Longitude 

-72.7449 

-73.1243 

-73.021 

-72.9881 

-73.0789 

-73.5908 

-74.1988 

Starting  Latitude 

31.99 

32.0734 

31.9384 

31.7897 

32.5184 

33.0707| 

33.887 

Ending  Longitude 

-72.7407 

-73.1082 

-73.019 

-73.0147 

-73.1451 

-73.7004 

-74.2319 

Ending  Latitude 

31.9744 

32.0842 

31.9196 

31.7416 

32.6475 

33.2104 

34.0078 

Figure  4  T-REMUS  deployment  times  and  tracks. 


From  14  June  to  18  June  our  3  buoy  system  drifted  into  an  increasingly  stronger  mesoscale  field  with 
both  shear  and  strain  variability  increasing.  Figure  5  shows  a  satellite  derived  map  of  SST  and  our 
Gateway  buoy  tracks. 


SST  [°  C]  14-Jun-201 1 
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Figure  5  Track  (dotted  circles)  of  the  3  buoy  system  from  14  June  2011  to  18  June  2011 

overlaid  on  SST  map. 
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In  figure  6  we  show  a  contour  map  of  salinity  obtained  from  the  T-REMUS  CTD  on  18  June  2011.  The 
vertical  axis  is  density  and  the  horizontal  axis  is  the  along  track  location  of  the  observations  in  a  fixed 
earth  coordintate  system.  Note  that  density  for  this  environment  is  mainly  determined  by  temperature 
so  that  salinity  can  be  used  approximately  as  a  passive  tracer,  although  there  are  some  instances  where 
double  diffusion  effects  can  occur.  Figure  6  clearly  shows  three  distinct  regions,  which  are  labelled  in 
the  Figure.  Region  2  indicates  a  cold  fresh  intrusion.  Toward  the  end  of  the  run  starting  at 
approximately  10  km  downstream  frontal  strucutre  is  observed. 


Salinity  (PSU)  18  June  2011 


downstream  distance  (km) 


Figure  6  contour  map  of  Salinity  as  a  function  of  density  and  downstream  distance  obtained  from 

the  t-REMUS  CTD. 


In  Figure  7  we  show  a  contour  map  of  %  the  rate  of  molecular  diffusion  of  temperature  (heat).  Note 
how  the  larges  values  are  typically  along  the  edges  of  the  intrusions, i.e  region  #  2. 
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Chi  IS  June  2011 


distance  (km) 


Figure  7  Contour  map  of  %  as  a  function  of  density  and  downstream  distance 


Following  the  approach  of  D'Asaro  (3008)  we  have  examined  the  along  isopycnal  mixing  using  a 
density  coordintate  system.  The  unqiue  combination  of  data  obtained  by  the  T-REMUS  vehilce,  e,x 
and  finescale  gradients,  both  along  and  across  isopycnals  of  T  ansd  S  allows  us  to  calculate  a 
temprature  budget  in  density  coordintates  ( McDougall,  1984)  This  is  shown  in  Figure  8.  We  see  that 
the  along  isopycnal  diffusivity  can  reach  vaues  of  order  10  m2/sec,  a  value  comparable  at  times  with 
the  with  LatMix  201 1  dye  estimate  during  the  same  time  of  deployment.  Thus  our  tentative  conclusion 
to  date  is  that  samll  scale  intrusions  (  vertical  sclaes  of  10-15  m  and  horizontal  scales  of  order 
kilometers)  can  play  a  very  signficant  role  in  along  isopycnal  mixing  and  should  be  taken  into  acount 
in  numerical  models. 
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Figure  8.  Top  the  advction  diffusion  equation  in  density  coordinates.  Bottom  3  paels  indicat:  (the 
dissipation  term  (shown  in  yellow  in  the  above  equation)  and  the  along  isopycnal  diffusion  term  (I- 
diff).  The  lower  two  panels  are  the  along  and  across  ispoycnal  diffusivities. 
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